Shell Corrections of Superheavy Nuclei in Self-Consistent Calculations 
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Shell corrections to the nuclear binding energy as a measure of shell effects in superheavy nuclei are 
studied within the self-consistent Skyrme-Hartree-Fock and Relativistic Mean-Field theories. Due 
to the presence of low-lying proton continuum resulting in a free particle gas, special attention is 
paid to the treatment of single-particle level density. To cure the pathological behavior of shell 
correction around the particle threshold, the method based on the Green's function approach has 
been adopted. It is demonstrated that for the vast majority of Skyrme interactions commonly 
employed in nuclear structure calculations, the strongest shell stabilization appears for Z=124, and 
126, and for A'^=184. On the other hand, in the relativistic approaches the strongest spherical 
shell effect appears systematically for Z=120 and A'^=172. This difference has probably its roots in 
the spin-orbit potential. We have also shown that, in contrast to shell corrections which are fairly 
independent on the force, macroscopic energies extracted from self-consistent calculations strongly 
depend on the actual force parametrisation used. That is, the A and Z dependence of mass surface 
when extrapolating to unknown superheavy nuclei is prone to significant theoretical uncertainties. 

PACS number(s): 21.10.Dr, 21.10.Pc, 21.60. Jz, 27.90. +b 



I. INTRODUCTION 

The stability of the heaviest and superheavy elements 
has been a long-standing fundamental question in nu- 
clear science. Theoretically, the mere existence of the 
heaviest elements with Z>10A is entirely due to quan- 
tal shell effects. Indeed, for these nuclei the shape of 
the classical nuclear droplet, governed by surface ten- 
sion and Coulomb repulsion, is unstable to surface dis- 
tortions driving these nuclei to spontaneous fission. That 
is, if the heaviest nuclei were governed by the classical 
liquid drop model, they would fission immediately from 
their ground states due to the large electric charge. How- 
ever, in the mid-sixties, with the invention of the shell- 
correction method, it was realized that long-lived super- 
heavy elements (SHE) with very large atomic numbers 
could exist due to the strong shell stabilization ||l|-y. 

In spite of tremendous experimental effort, after about 
thirty years of the quest for superheavy elements, the 
borders of the upper-right end of the nuclear chart are 
still unknown [pj. However, it has to be emphasized 
that the recent years also brought significant progress 
in the production of the heaviest nuclei pjq]. Dur- 
ing 1995-96, three new elements, Z=110, 111, and 112, 
were synthesized by means of both cold and hot fu- 
sion reactions |7|-|lC|]. These heaviest isotopes decay pre- 
dominantly by groups of a particles (a chains) as ex- 
pected theoretically ||ll|-[l3|. Recently, two stunning dis- 
coveries have been made. Firstly, hot fusion experi- 



ments performed in Dubna employing ^^Ca-f^'^^Pu and 
^^Ca-|-^'*^Pu "hot fusion" reactions ||3 gave evidence for 
the synthesis of two isotopes (^=287 and 289) of the 
element Z=114. Secondly, the Berkeley-Oregon team, 
utihzing the "cold fusion" reaction ^"^Kr-h^^^Pb [15), ob- 
served three a-decay chains attributed to the decay of 
the new element Z=118, A— 293. The measured a-decay 
chains ^^^114 and ^^^118 turned out to be consistent 
with predictions of the Skyrme-Hartree-Fock (SHF) the- 
ory ||l^ and the Relativistic Mean-Field (RMF) theory 

The goal of the present work is to study shell closures 
in SHE. To that end we use as a tool microscopic shell 
corrections extracted from self-consistent calculations. 
For medium-mass and heavy nuclei, self-consistent mean- 
field theory is a very useful starting point |lq| . Nowadays, 
SHF and RMF calculations with realistic effective forces 
are able to describe global nuclear properties with an 
accuracy which is comparable to that obtained in more 
phenomenological macroscopic-microscopic models based 
on the shell-correction method. 

In previous work [nol, shell energies for SHE elements 
were extracted by subtracting from calculated HF bind- 
ing energies the macroscopic Yukawa-plus-exponential 
mass formula [Q with parameters of Ref. |2^. In an- 
other work, based on the RMF theory pi], shell cor- 
rections were extracted for the heaviest deformed nu- 
clei using the standard Strutinsky method in which the 
positive-energy spectrum was approximated by quasi- 
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FIG. 1. Single-proton levels for A'^=184 isotones with 110<Z<130 calculated in the Skyrme-Hartree-Fock model with SkM* 
(left) and SkP effective interactions. Positive (negative) parity levels are indicated by solid (dashed) lines and by their spherical 
labels (nlj). Note that in both cases the nucleus Z=126 is proton-unbound, i.e., the pi/2 shell has positive energy. 



bound states. Neither procedure can be considered as 
satisfactory. A proper treatment of continuum states is 
achieved with a Green's function method |2^. Wc cm- 
ploy this method for the present study of shell corrections 
of SHE. 

The material contained in this study is organized as fol- 
lows. T he m otivation of this work is outlined in Sect. O, 
Section III contains a brief discussion of the Strutinsky 
energy theorem on which the concept of shell correc- 
tion is based. The Green's Function HP method used 
to extract the single-particle level density is presented 
in Sect. IV. Section discusses the details of our HF 
and RMF models and describes the Strutinsky proce- 
dure employed. The results of calculations for shell cor- 
rections in spherical SHE and for macroscopic energies 
extracted from self-consistent binding energies are dis- 



state, is in nice agreement 
p|plp^,|2g|. Still heavier 



cussed in Sec. VI. Finally, Sec. VII contains the main 



conclusions of this work. 



II. MOTIVATION 

All the heaviest elements found recently are believed 
to be well deformed. Indeed, the measured a-decay 
energies, along with complementary syntheses of new 
neutron-rich isotopes of elements Z=106 and Z=108, 
have furnished confirmation of the special stability of 
the deformed shell at N—162 predicted by theory P4| , p5| . 
Beautiful experimental confirmation of large quadrupole 
deformations in this mass region comes from gamma- 
ray spectroscopy. Recent experimental works P6lp7| suc- 
ceeded in identifying the ground-state band of ^^No (the 
heaviest nucleus studied in gamma-ray spectroscopy so 
far). The quadrupole deformation of ^'^'*No, inferred from 



the energy of the deduced 2^^ 
with theoretical predictions 
and more neutron-rich elements are expected to be spher- 
ical due to the proximity of the neutron shell at A^=184. 
This is the region of SHE which we will investigate here. 
In spite of an impressive agreement with available ex- 
perimental data for the heaviest elements, theoretical un- 
certainties are large when extrapolating to unknown nu- 
clei with greater atomic numbers. As discussed in Refs. 
||l9| , p0| , the main factors that influence the single-proton 
shell structure of SHE are (i) the Coulomb potential and 
(ii) the spin-orbit splitting. As far as the protons are 
concerned, the important spherical shells are the closely 
spaced lii3/2 and 2/7/2 levels which appear just below 
the Z—114 gap, the 2/5/2 shell which becomes occupied 
at Z—120, the 3p3/2 shell which becomes occupied at 
Z—124, and the 3pi/2 and liii/2 orbitals whose split- 
ting determines the size of the Z—126 magic gap. In- 
terestingly, while the ordering of single-proton states is 
practically the same for all the self-consistent approaches 
with realistic effective interactions (see Fig. |l|and single- 
particle diagrams in Refs. Il9,pfl]), their relative posi- 
tions vary depending on the choice of force parameters. 
Since in the region of SHE the single-particle level den- 
sity is relatively large, small shifts in positions of single- 
particle levels can influence the strength of single-particle 
gaps and be crucial for determining the shell stability of 
a nucleus. As a result, there is no consensus between 
theorists concerning the next proton magic gap beyond 
Z—82. While most macroscopic-microscopic (non-self- 
consistent) approaches predict Z~ll4 to be magic, self- 
consistent calculations suggest that the center of the pro- 
ton shell stability should be moved up to higher proton 
numbers, Z=120, 124 or 126 H&lsl. It is to be noted 
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FIG. 2. Single-neutron levels for Z=120 isotopes with 160<Af<190 calculated in the RMF approach with NL3 (left) and 
NL-Z2 non-linear parametrizations. The line convention is the same as in Fig. PM. Note large neutron gaps at N—172 and 184. 



that the Coulomb potential mainly influences the magni- 
tude of the Z=11A gap. (Here, the self-consistent treat- 
ment of Coulomb energy is a key factor.) On the other 
hand, the spin-orbit interaction determines the position 
of the 2/ and ip shells which define the proton shell struc- 
ture above Z>114. 

The spherical neutron shell structure is governed by 
the following orbitals: IJ15/2 (below the N—IQA gap), 
257/2, 3^5/2, 3(^3/2, and 4si/2 and IJ13/2 whose splitting 



determines the size of the N 
and Refs. ^^). 



184 spherical gap (see Fig. 
Again, similar to the proton case, the 
order of the single-neutron orbitals between A^=164 and 
184 is rather robust, while sizes of single-particle gaps 
vary. For instance, the iV=172 gap, predicted by the 
RMF calculations shown in Fig. I2I results from the large 
energy splitting between the 257/2 and 8^5/2 shells. In 
non-relativistic models, these two orbitals are very close 
in energy, and this degeneracy is related to the pseudo- 
spin symmetry [ p2|j3^ . Interestingly, in the SHF calcu- 
lations, the pseudo-spin degeneracy holds in most cases. 
Namely, certain neutron orbitals group in pairs (pseudo- 



(3d; 



'Z/2; 



As 



I/2J 



and the 



spin doublets): (257/2, 3^5/2^ 

same holds for proton orbitals, e.g., (2/5/2, 3^3/2). Con- 
sidering the fact that the idea of pseudo-spin has rela- 
tivistic roots fsj,^ , it is surprising to see that this sym- 
metry is so dramatically violated in the RMF theory. As 
a matter of fact, the presence of pronounced magic gaps 
at Z=120 and iV=172 in RMF models (see below) is a 
direct manifestation of the pseudo-spin symmetry break- 
ing. 

As discussed in Ref. [19l , neutron-deficient superheavy 
nuclei are expected to be unstable to proton emission. In- 
deed, as seen in Fig. || the proton 3pi/2 shell has positive 
energy for Z>126, i.e., in these nuclei the 3pi/2 level is 
a narrow resonance. Due to huge Coulomb barriers, su- 



perheavy nuclei with Qp<1.5 MeV are practically proton- 
stable |19]. However, the higher-lying single-proton or- 
bitals are expected to have sizable proton widths. 

In order to assess the magnitude of shell effects de- 
termined by the bunchiness of single-particle levels, it is 
useful to apply the Strutinsky renormalization procedure 
p6|-p8| which makes it possible to calculate the shell cor- 
rection energy. Unfortunately, the standard way of ex- 
tracting shell correction breaks down for weakly bound 
nuclei where the contribution from the particle contin- 
uum becomes important p9| . Recently, a new method of 
calculating shell correction, based on the correct treat- 
ment of resonances, has been developed |^,^. The im- 
proved method is based on the theory of Gamow states 
(eigenstates of one-body Hamiltonian with purely outgo- 
ing boundary conditions) which can be calculated numer- 
ically for commonly used optical- model potentials p2[ . 
While this "exact" procedure cannot be easily adopted 
to the case of microscopic self-consistent potentials, its 
simplified version applying the Green's-function method 
can H. 



III. SHELL CORRECTION AND THE ENERGY 
THEOREM 

The main assumption of the shell-correction (macrosco- 
pic-microscopic) method [|6-38 4^ is that the total en- 
ergy of a nucleus can be decomposed into two parts: 



E = E + E. 



shell; 



(1) 



where E is the macroscopic energy (smoothly depend- 
ing on the number of nucleons and thus associated with 
the "uniform" distribution of single-particle orbitals) and 



-E'shcU is the shell-correction term that fluctuates with 
particle number reflecting the non-uniformities (bunchi- 
ness) of the single-particle level distribution. In order to 
make a separation (|^), one starts from the one-body HF 
density matrix p 



p{r', 






nMr')^:{r) 



(2) 



which can be decomposed into a "smoothed" density p 
and a correction 5p^ which fluctuates with the shell filling 



p = p + Sp. 



(3) 



In Eq. (y), Tii is the single-particle occupation coefficient 
which is equal to 1(0) if the level e^ is occupied (empty). 
The smoothed single-particle density p can be expressed 
by means of the smoothed distribution numbers hi [H : 



Kr',r) ^y^jL,<j),{r')(t)*{r). 



(4) 



When considered as a function of the single-particle en- 
ergies Ci, the numbers hi vary smoothly in an energy in- 
terval of the order of the energy difference between major 
shells. The averaged HF Hamiltonian hup can be directly 
obtained from p. The expectation value of a HF Hamil- 
tonian (containing the kinetic energy, t and the two-body 
interaction, v) can then be written in terms of p and 6p 

MM- 

Ehf - Tr(tp) + ^Tr Tiipvp) =E + E^sc + 0{Sp^), (5) 



where 



E — Tr(t/5) + -Tr Tiipvp) 
is the average part of E'hf and 



(6) 



Eosc = Tr(/iHF^p) with /ihf =t + Tr(i;p) (7) 

is the first-order term in dp representing the shell- 
correction contribution to E^p. If a deformed phe- 
nomenological potential gives a similar spectrum to the 
averaged HF potential /ihf, then the oscillatory part of 
-Ehf, given by Eq. (0), is very close to that of the de- 
formed shell model, £'shcii=-£'osc + 0((5p^). The second- 
order term in Eq. (H) is usually very small and can be 
neglected [^. The above relation, known as the Struti- 
nsky Energy Theorem, makes it possible to calculate 
the total energy using the non-self-consistent, deformed 
independent-particle model; the average part E is usually 
replaced by the corresponding phenomenological liquid- 
drop (or droplet) model value, -Emacro- It is important 
that -Eshcii must not contain any regular (smooth) terms 
analogous to those already included in the phenomeno- 
logical macroscopic part. The numerical proof of the En- 
ergy Theorem was carried out by Brack and Quentin p7| ] 
who demonstrated that Eq. (|^) holds for i^shcU defined 



by means of the smoothed single-particle energies (eigen- 
values of /ihf)- 

In this work, we use a simpler expression to extract 
the shell correction from the HF binding energy, which 
should also be accurate up to 0((5/9^). Namely, as 
an input to the Strutinsky procedure we take the self- 
consistent single-particle HF energies, ef^. In this case, 
the shell correction is given by 



Eshcwip) = '^{^i - hi)ei + 0{5p^). 



(8) 



The equivalent macroscopic energy can easily be com- 
puted by taking the difference 



-Emacro ~ E = E{p) — £'shcll(p)- 



(9) 



IV. GREEN'S FUNCTION HARTREE-FOCK 
APPROACH TO SHELL CORRECTION 

The HF equation is generally solved using a harmonic 
oscillator expansion method or by means of a discretiza- 
tion in a three-dimensional box. In both cases, a great 
number of unphysical states with positive energy appear. 
The effect of these quasi-bound states is disastrous for 
the Strutinsky renormalization procedure p^ , p9| - ^ , p8[ . 
Indeed, if one smoothes out the single-particle energy 
density. 



ffsp 



(e)-E'5' 



„HF 



). 



(10) 



it would diverge at zero energy because the presence of 
the unphysical positive energy states. Consequently, the 
resulting shell correction becomes unreliable. 

In order to avoid the divergence of g{e) around 
the threshold, we apply the Green's-function method 
p3| , [49| -|5^ for the calculation of the single-particle level 
density. In this method, the level density is given by the 
expression 



lie) 



i3{Tr[G+(e)-G+,,(e)]} 



(11) 



where G^{e) — (e — h + iO) ^ is the outgoing Green's op- 
erator of the single-particle Hannltonian h{p), and Gf^.^^ 
is the free outgoing Green's operator that belongs to the 
"free" single-particle Hamiltonian. This latter is derived 
from the full HF Hamiltonian in such a way that those 
terms are kept which are related to the kinetic energy 
density and to the direct Coulomb term. The interpre- 
tation of Eq. (ni]) is straightforward: the second term in 
Eq. (11) contains the contribution to the single-particle 
level density originating from the gas of free particles. 

The single-particle level density defined by the Green's- 
function expression (O) behaves smoothly around the 
zero-energy threshold; for finite-depth Hamiltonians this 



definition is the only meaningful way of introducing 5(e). 
The level density (|l^) automatically takes into account 
the effect of the particle continuum which r nav influence 
the results of shell-correction calculations [ [JOplJI , espe- 
cially pronounced for systems where the Fermi level is 
close to zero, i.e., drip- line nuclei. 

Because it is difficult to calculate the Green's-function, 
in this work we applied the approximation introduced 
in Ref. [E3[. In this approach, the single-particle level 
density is expressed as 



ff(e) 



,HF' 






-E^ 



HF.froo 



(12) 



where e^ ' '^^^ are the eigenvalues of the free one-body 
HF Haniiltonian. As usual in the Strutinsky procedure, 
the smooth level density can be obtained by folding g(e) 
with a smoothing function f{x): 



1 



/ <J —00 



+ 00 



de' g{e') f 



5o(e) -5frco(e), 



(13) 



where 7 is the smoothing width, go{e) is the smooth level 
density obtained from the HF spectrum (including the 
quasi-bound states), and girocie) is the contribution to 
the smooth level density from the particle gas. 

In practice, g{e) can be calculated in three steps. 
First, we solve the HF equations to determine the self- 



consistent energies e 



HF 



In the next step, we calcu- 



late the positive-energy gas spectrum e^ ' ^'^'^ at the self- 
consistent minimum. In particular, we take the Coulomb 
force from the self-consistent calculation. Finally, we 
compute go(e) and ^free(e) using the same folding func- 
tion. The quality of approximation (n2|) was tested 
in Ref. [E3[ where it was demonstrated that, when in- 
creasing the number of basis states, the resulting single- 
particle level density quickly converges to the exact re- 
sult. 



V. SELF-CONSISTENT MODELS 

A. Skyrme-Hartree-Fock Model 

In the SHF method, nucleons are described as nonrela- 
tivistic particles moving independently in a common self- 
consistent field. Our implementation of the HF model is 
based on the standard ansatz [Q . The total binding en- 
ergy of a nucleus is obtained self-consistently from the 
energy functional: 



£ = £ki 



kin 



£. 



Sk 



£c 



f £sk,ls 
f pair ^ f CM , 



(14) 



where £^kin is the kinetic energy functional, £sk is the 
Skyrme functional, £sk,is is the spin-orbit functional, £c 



is the Coulomb energy (including the exchange term), 
^pair is the pairing energy, and f cm is the center-of-mass 
correction. 

Since there are more than 80 different Skyrme param- 
eterizations on the market, the question arises, which 
forces should actually be used when making predictions 
and comparing with the data? Here, we have chosen 
a small subset of Skyrme forces which perform well for 
the basic ground-state properties (masses, radii, sur- 
face thicknesses) and have sufficiently different proper- 
ties which allows one to explore the possible variations 
among parameterizations. This subset contains: SkM* 
g, SkT6 H], Z„ g, SkP I?), SLy4 |§, and Skll, 
SkI3, and SkM from Ref. |9|. We have also added the 
force SkO from a recent exploration [Q . Most of these 
interactions have been used for the investigation of the 
ground-state properties of SHE before |p|,^j29|-|3l|] . All 
the selected forces perform well concerning the total en- 
ergy and radii. They all have comparable incompressibity 
i4r=210-250MeV and comparable surface energy which 
results from a careful fit to ground-state properties Q . 
Variations occur for properties which are not fixed pre- 
cisely by ground-state characteristics. The effective nu- 
cleon mass is 1 for SkT6 and SkP, 0.9 for SkO, around 0.8 
for SkM* and Z„, and even lower, around 0.65, for SLy4, 
Skll, Skis, and SkM. Isovector properties also exhibit 
large variations. For SkI3 and SkM, the spin-orbit func- 
tional is given in the extended form of [p9| which allows a 
separate adjustment of isoscalar and isovector spin-orbit 
force. The standard Skyrme forces use the particular 
combination of isoscalar and isovector terms which were 
motivated by the derivation from a two-body zero-range 
spin-orbit interaction |6^]. (For a detailed discussion of 
the spin-orbit interaction in SHF we refer the reader to 
Refs. |0|J59|j6|-||.) 



B. Relativistic Mean-Field Model 

In our implementation of the RMF model, nucleons are 
described as independent Dirac particles moving in lo- 
cal isoscalar-scalar, isoscalar- vector, and isovector-vector 
mean fields usually associated with a, lu, and p mesons, 
respectively pq] . These couple to the corresponding local 
densities of the nucleons which are bilinear covariants of 
the Dirac spinors similar to the single-particle density of 
Eq. (I). 

The RMF is usually formulated in terms of a covariant 
Lagrangian; see, e.g., Ref. |^. For our purpose we prefer 
a formulation in terms of an energy functional that is 
obtained by eliminating the mesonic degrees of freedom 
in the Lagrangian. For a detailed discussion of the RMF 
as an energy density functional theory, see Refs. [|66|-|69|. 
The energy functional of the nucleus 



'RMF 



-— £]<i 



kin 



£a 

£c 



^UJ t Oq 



-CM 



(15) 



is composed of the kinetic energy of the nucleons £^kin, 
the interaction energies of the a^ uj^ and p fields, and 
the Coulomb energy of the protons Eq. All these are 
bilinear in the nucleonic densities as in the case of non- 
relativistic models [cf. Eq. (H)]. Pairing correlations arc 
treated in the BCS approach employing the same non- 
relativistic pairing energy functional f pair that is used in 
the SHF model. The center-of-mass correction £cm is 
also calculated in a non-relativistic approximation; see 
|70| for a detailed discussion. The single-particle ener- 
gies Ci needed to calculate the shell correction are the 
eigenvalues of the one-body Hamiltonian of the nucleons 
which is obtained by variation of the energy functional 

®- 
In the context of our study, it is important to note that 

the spin-orbit interaction emerges naturally in the RMF 
from the interplay of scalar and vector fields 1 65 1 . With- 
out any free parameters fitted to single-particle data, the 
RMF gives a rather good description of spin-orbit split- 
tings throughout the chart of nuclei |^ . 

As in SHF, there exist many RMF parameterizations 
which differ in details. For the purpose of the present 
study, we choose the most successful (or most commonly 
used) ones: NLl 0, NL-Z |7|], NL-Z2 [|o), NL-SH @, 
NL3 Ifri, and TMl |75l. AU of them have been used for 

11,0 • 



investigations of SHE 




The parameterization NLl is a fit of the RMF along 
the strategy of Ref. ||5^ used also for the Skyrmc inter- 
action Zcr. The NL-Z parametrization is a refit of NLl 
where the correction for spurious center-of-mass motion 
is calculated from the actual many-body wave function, 
while NL-Z2 is a recent variant of NL-Z with an improved 
isospin dependence. The force NL3 stems from a fit in- 
cluding exotic nuclei, neutron radii, and information on 
giant resonances. The NL-SH parametrization was fitted 
with a bias toward isotopic trends and it also uses infor- 
mation on neutron radii. The force TMl was optimized 
in the same way as NL-SH except for introducing an addi- 
tional quartic self-interaction of the isoscalar- vector field 
to avoid instabilities of the standard model which occur 
for small nuclei. For SHE, the results obtained with NL- 
Z are not distinguishable from results obtained with the 
parameterization PL-40 that is contained in exactly the 
same manner as NL-Z but uses a stabilized non-linearity 
of the scalar- isoscalar field [Q . (PL-40 was employed in 
some recent investigations of the properties of superheavy 
nuclei ||j3l|,0.) 

All the above parameterizations provide a good de- 
scription of binding energies, charge radii, and surface 
thicknesses of stable spherical nuclei with the same over- 
all quality as the SHF model. The nuclear matter prop- 
erties of the RMF forces, however, show some system- 
atic differences as compared to Skyrme forces. All RMF 
forces have comparable small effective masses around 
m* /m K 0.6. (Note that the effective mass in RMF 
depends on momentum; hence the effective mass at 
the Fermi energy is approximately 10% larger.) Com- 
pared with the SHF model, the absolute value of the 



energy per nucleon is systematically larger, with val- 
ues around —16.3 MeV, while the saturation density is 
always slightly smaller with typical values around 0.15 
nucleons/fm'^. The compressibility of the RMF forces 
ranges from low values around 170 MeV for NL-Z to 
K~?>^^ MeV for NL-SH, which is rather high. There are 
also differences in isovector properties; the symmetry en- 
ergy coefficient of all RMF forces is systematically larger 
than for SHF interactions, with values between 36.1 MeV 
for NL-SH and 43.5 MeV for NLl (see discussion below). 



C. Details of Calculations 

In order to probe the single-particle shell structure of 
SHE, SHF and RMF calculations were carried out under 
the assumption of spherical geometry. By doing so we 
intentionally disregard deformation effects which make it 
difficult to compare different models and parametriza- 
tions. For the same reason, pairing correlations were 
practically neglected. (In order to obtain self-consistent 
spherical solutions for open-shell nuclei, small constant 
pairing gaps, A<100keV were assumed; the correspond- 
ing pairing energies are negligible. This procedure is ap- 
proximately equivalent to the filling approximation.) 

The SHF calculations were carried using the 
coordinate-space Hartree-Fock code of Ref. [|o|. The 
HE equations were solved by the discretization method. 
To obtain a proper description of quasi-bound states, 
it was necessary to take a very large box and a very 
dense mesh. The actual box size was chosen to be 21fm 
and the mesh spacing was 0.3 fm. With this choice, the 
low-lying positive-energy proton states obtained in SHF 
perfectly reproduce proton resonances obtained by solv- 
ing the Schrodinger equation for the HF potential with 
purely outgoing boundary conditions. 

The Strutinsky procedure contains two free parame- 
ters, the smoothing parameter 7 and the order of the cur- 
vature correction p. In calculating the Strutinsky smooth 
energy, instead of the traditional plateau condition we 
applied the generalized plateau condition described in 
Ref. [Bl|. The optimal values of 7 (in units of oscilla- 
tor frequency hjjQ—41/A^/^) calculated for several nuclei 
turned out to be close to 7p=1.54 and 7„=1.66 for pro- 
tons and neutrons, respectively; these values, together 
with p=10, were adopted in our calculations of shell cor- 
rections in SHF. 

In the RMF approach, the shell correction can be ex- 
tracted from the single-particle spectrum like in SHF. To 
demonstrate it, one proceeds along the steps discussed 
in Sec. [II. The total RMF energy ( p^ can be decom- 
posed into a smooth part and a correction that fluctuates 
according to the actual level density. Since the RMF en- 
ergy functional is bilinear in the densities, the extracted 
shell correction should be accurate up to order 0{5p^). 

The RMF calculations were carried out using the 
coordinate-space code of Ref. ISlI] . As in the SHF case. 
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FIG. 3. Spherical neutron shell corrections for the Z—120 isotopes calculated in nine Skyrme-Hartree-Fock models. The 
dotted line for SkM* shows the proton shell correction for comparison. In all nine cases, the minimum of the shell correction 
is predicted at A''=184. 



the box size was chosen to be 21 fm with a mesh spacing 
of 0.3fm. 

As already mentioned, all successful RMF parameter- 
izations give a rather small effective mass. This leads to 
a small level density around the Fermi surface which in 
turn requires a very large smoothing range 7 when cal- 
culating the smoothed level density g. The values for 7 
are strongly correlated with the order of the curvature 
correction polynomial p |41[| ; the value p=10 chosen here 
is large enough to provide in nearly all cases a sufficiently 
smooth g, but also small enough that we can restrict the 
model space to levels up to 60 MeV, which is much larger 
than the space used in usual RMF calculations. We have 
adjusted the smoothing range 7 to the actual level den- 
sity of a large number of nuclei to fulfill a generalized 
plateau condition along the strategy of pi. This leads 
always to values around 7p=2.0 for protons and 7„=2.2 
for neutrons. All results presented in this paper are cal- 
culated with p=10 and 7 fixed at these values. 



VI. RESULTS 

A. Spherical Shell Corrections in Superheavy Nuclei 

According to the SHF calculations of Rcf. |Q , the spher- 
ical magic neutron number in the SHE region is A'^=184; 
all the N—184: isotones have been predicted to have 
spherical shapes. The magicity of iV=184 in SHF is con- 
firmed in this study. Figure ^ displays neutron shell cor- 
rection calculated in several SHF models as a function 
of N for Z=120. The absolute minimum of shell energy 
always appears at A^=184. The A^=172 shell effect is 
also seen, but it exhibits a strong force-dependence (it is 
particularly pronounced for Z^, Skl3, SkM and SLy4). 

As already mentioned, the neutron levels have the 
same ordering for nearly all forces; all differences seen 
in the shell corrections are therefore caused by slight 
changes in the relative distances of the single-particle lev- 
els between the models. Forces with large effective masses 
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FIG. 4. Spherical proton shell corrections for the chains of Ai'=184 isotones (solid lines) and N=172 isotones (dotted lines) 
calculated in nine Skyrme-Hartree-Fock models. The dashed (dash-dotted) line for SkM* shows the neutron shell correction 
for A'^=184 {N=172). For A'^=184, the minimum of the shell correction is predicted at Z=124-126 for all parametrizations. 



like SkO, SkP, and SkT6, give a comparatively large level 
density which washes out the shell effects below A^=184. 
Forces with small effective masses (i.e., smaller level den- 
sity) are much more likely to show significant shell effects 
at lower neutron numbers around N—172. 

At fixed Z, the proton shell correction changes rather 
gradually as a function of neutron number; this is illus- 
trated in Fig. H for the Skyrme force SkM* . (Most of the 
Skyrme forces give a similar result.) Note that the proton 
shell corrections are generally smaller than those for the 
neutrons. At a second glance, however, one sees that the 
slow variations of the proton shell correction with neu- 
tron number are correlated with neutron shell closures. 
For instance, the Z=120 shell correction is largest at neu- 
tron numbers around N=172 and it becomes reduced 
when approaching A^=184. This is caused by the self- 
consistent rearrangement of single-particle levels accord- 
ing to the actual density distribution in the nucleus and 
cannot appear in macroscopic-microscopic models with 



assumed average potentials (see Refs. [ pOpl| ] for more 
discussion related to this point). 

Proton shell corrections for the A^=184 and iV=172 
isotones, obtained in the SHF model, are displayed in 
Fig. Q as a function of Z. For SkM* , neutron shell correc- 
tions are also shown for the 7V=172 and iV=184 isotones. 
The shift of the magic proton number with neutron num- 
ber when going from A^=172 to iV=184 is clearly visi- 
ble. For A^=172 most of the Skyrme forces (exceptions 
are SkT6 and SkP) agree on a magic Z — 120, while 
for N = 184 the shell correction shows a minimum at 
Z= 124-126 in all cases. (Actually, in most cases, shell- 
corrections slightly favor Z=124 over Z—126; this is re- 
lated to the gradual increase of single-particle energies of 
3p3/2 and 3pi/2 orbitals above Z—120.) 

Proton shell corrections and the iV=172 neutron shell 
corrections are systematically smaller than those for neu- 
trons at A''=184. This partly explains why spherical 
ground states of SHE are so well correlated with the 
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FIG. 5. Spherical neutron shell corrections for the Z—12Q isotopes calculated in six relativistic mean-field models. In all six 
cases, the minimum of the shell correction is predicted at N—172. For NL-Z2 the proton shell correction is given for comparison 
by the dashed line. 



magic neutron number A''=184, see, e.g., p^ , p2| , p9[ . Note 
that for the majority of Skyrme forces the N=172 iso- 
tones are predicted to be deformed. 

Skyrme forces with non-standard isospin dependence 
of the spin-orbit interaction are the only ones that give 
additional (but not very pronounced) shell closures. In 
the Skl4 model, there appears a secondary minimum 
at Z==114 for A^=184, while Skl3 is the only Skyrme 
force which points at Z—120 also for A^=184. A non- 
standard spin-orbit interaction, however, does not necce- 
sarily lead to shell closures other than Z=124-126 for 
N=18A. For SkO, which has a spin-orbit force that is 
similar to SkI4, the Z=114 shell is only hinted. It is to 
be noted that for several interactions such as Z^, Sklx, 
and SkO, shell correction changes rather slowly between 
Z—llA and Z=126. This indicates that none of the pro- 
ton shell gaps in this region can be considered as truly 
"magic" . (The weak Z-dependence of proton shell correc- 
tion above Z—\1A was pointed out in the early Ref. [p2|.) 

The RMF results presented in Figs. H and show a 
pattern that is internally consistent but different from 
that of SHF. The minimum of neutron shell correction is 
systematically predicted at N=n2. Except for NL-SH 
and TMl, the shell effect at 7V=182-184 is also clearly 
seen. Note that the 7V=184 gap in the single-particle 
spectrum is in all cases larger than the one at A^=182 
(see Fig. ||). The gaps are separated by a single 4si/2 
level which contributes very weakly to the shell energy. 
To illustrate the variation of proton shell effects along 



the Z=120 chain, proton shell corrections in NL-Z2 are 
also displayed in Fig. ^. Their pattern is very similar to 
that obtained in SHF models. 

Looking at the proton shell corrections along the chain 
of A^=184 isotones, see Fig. 0, the strongest shell effect 
is now obtained for Z=120. When comparing the results 
for the A^=184 and A^=172 chains, it can be seen again 
that the proton shell correction at Z=12Q is strongly 
correlated with neutron number iV=172. However, un- 
like in SHF, the ^=120 shell does not vanish completely 
for N—184:. Proton shell corrections obtained with NLl, 
NL-Z, and NL-Z2 at N=18A vary rather slowly between 
Z—120 and Z=126, and this resembles the patterm ob- 
tained in SHF. Again, as in the case of Skyrme forces, 
proton shell corrections in RMF are smaller than those 
for the neutrons (cf. NL-Z2 calculations in Fig. o). The 
increase in the proton shell correction at very large val- 
ues of Z for TMl is related to the spherical Z=132 shell 
predicted by this interaction |pi|| . 

Shell closures can also be analysed in terms of the two- 
neutron and two-proton shell gaps 

S2n - E{N + 2,Z)- 2E{N, Z) + E[N ~ 2, Z), 

62p - E{N, Z + 2)- 2E{N, Z) -f E{N, Z-2) (16) 

discussed in Re fs. [^ 77|. The pattern of shell corrections 
calculated in SHF and RMF nicely follows the behavior of 
neutron and proton shell gaps found there. In particular, 
the strong correlation between shell effects at Z=12Q and 
N~172 in RMF is seen in both representations. While 
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FIG. 6. Spherical proton shell corrections for the Ai'=184 (solid line) and N=172 (dotted line) isotones calculated in six 
relativistic mean-field models. For all parametrizations, the minimum of the proton shell correction is predicted at Z=120. 
For NL-Z2 the neutron shell correction for the A'^=184 (dashed line) and N=172 (dash-dotted line) isotones are given for 
comparison. 



shell gaps are related (but not equivalent) to the gaps 
in the single-particle spectrum, the shell correction gives 
also a measure of the stabilizing effect of a shell closure 
on the nuclear binding energy. 



B. Macroscopic Energies 

By subtracting the shell correction from the calculated 
binding energy, one obtains a rough estimate for the asso- 
ciated macroscopic energy -Emacro @- The macroscopic 
part of the SHF and RMF energies for the A^=184 iso- 
tones as a function of Z is displayed in Fig. ^ The macro- 
scopic energy of the Yukawa-plus-exponential mass for- 
mula of the Finite- Range Liquid Drop Model (FRLDM) 
of Rcf . 1221 , with parameters of Ref . |2^ , is also shown 
for comparison. To illustrate .Z-dependence, all energies 
were normalized to the value at Z=100. In general, the 
behavior of i^macro is similar in all cases. In particular, 
the macroscopic proton drip line is consistently predicted 
to be at ZRil20-124. It is interesting to note that the only 
Skyrme force which agrees with FRLDM is SLy4; other 
forces deviate from it significantly. The RMF forces give 
qualitatively the same results; there are several forces 
(NL-Z, TMl, and NL-SH) which give values of E^e.cro 
close to the FRLDM. 

In an attempt to understand the pattern shown in 
Fig. R, we employed the simple liquid drop model ex- 



pression 

-Emacro,LDM = Ovol^ + ^surf ^ 

{N - Zf 

-r "sym , 



(17) 



acoui 



yll/3- 



The parameters ai of Skyrme and RMF forces were cal- 
culated in the limit of symmetric nuclear matter; they are 
given in Table |, together with the values for the stan- 
dard liquid drop model (LDM) of Ref. ^ . [Note that 
these values slightly change when including higher-order 
terms in the LDM expansion (|l7|).] Figure |g shows the 
macroscopic energy (O) as a function of Z for the A^=184 
isotones. The huge differences between results for vari- 
ous Skyrme and RMF parametrizations can be traced 
back to their different symmetry-energy coefficients. In- 
deed, for most of the forces discussed, Osym is signifi- 
cantly greater than that of LDM, and this results in an 
increased slope of ii'macro.LDM. For the RMF forces the 
significantly larger avoi even further increases the differ- 
ence with respect to the LDM. Unfortunately, there is 
very little similarity between the results of microscopic 
calculations of Fig. and the results of expansion (|lj). 
When comparing the energy scales of Figs. and 0. one 
finds huge differences, of the order of 100 MeV, between 
£^macro and ^^macroXDM- While for RMF the energy or- 
dering remains the same in both cases, this feature does 
not hold for SHF. Only when looking at -Emacro.LDM, the 
results are ordered according to the corresponding val- 
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FIG. 7. Macroscopic energy E (91) extracted from the calculated Hartree-Fock energies of the 7V=184 isotones. For 
comparison, the phenomenological macroscopic energy of the Yukawa-plus-exponential mass formula (FRLDM) of Ref. |2d] 
with parameters of Ref. pl| is also displayed. To illustrate Z-dependence, all energies were normalized to the FRLDM value 
at Z=100. 



ues of asym, as expected. All this indicates that even for 
very heavy nuclei with A~300, the simple leptodermous 
expansion with parameters taken from nuclear matter 
calculations is not going to work ^,Q; the finite-size 
effects are still very important for SHE. 

In spite of the fact that macroscopic energies extracted 
from different self-consistent models systematically dif- 
fer, the corresponding shell corrections are similar. For 



TABLE L Key properties of symmetric nuclear matter for 
the Skyrme and RMF forces used in this paper: binding en- 
ergy per nucleon, surface energy, and symmetry energy, all 
in MeV. The RMF values for a^urf are taken from Ref. |||. 
The standard liquid drop model (LDM) values |84] are also 
shown. 



Force 


fflvol 


Hsurt 


Q-sym 


SkM* 


-15.9 


17.59 


30.0 


z. 


-15.9 


16.94 


26.7 


SkT6 


-16.1 


18.12 


29.9 


SLy4 


-16.1 


18.18 


32.0 


Skll 


-15.9 


17.31 


37.5 


SkI3 


-16.0 


17.52 


34.8 


SkM 


-15.9 


17.28 


29.5 


SkP 


-16.0 


17.95 


30.0 


SkO 


-15.8 


17.00 


32.0 


NLl 


-16.4 


18.66 


43.5 


NL-Z 


-16.2 


17.72 


41.7 


NL-Z2 


-16.1 




39.0 


NL3 


-16.2 


18.46 


37.4 


NL-SH 


-16.3 


19.05 


36.1 


TMl 


-16.3 




36.9 


LDM 


-15.7 


18.56 


28.1 



instance, the general pattern and magnitude of shell en- 
ergies displayed in Figs. || and ^ do not depend very much 
on the Skyrme interaction used, and the same is true for 
the RMF results shown in Figs, g and g. This means that 
although the global properties of effective interactions 
employed in this work differ, their single-particle spec- 
tra are fairly similar. Hence, shell corrections extracted 
from self-consistent single- particle spectra are very use- 
ful measures of spectral properties of effective forces. 
Figure [^ also illustrates how dangerous it is to extrap- 
olate self-consistent results in the region of SHE. The 
trends of relative binding energies (e.g., Qa values) are 
expected to smoothly deviate from force to force. The 
nice agreement with experimental data for the heaviest 
elements obtained in the SHE calculations with SLy4 |lq ] 
and in the macroscopic-microscopic calculations with the 
FRLDM 1 21 1 indicates that the macroscopic energies of 
forces which are too far off the FRLDM values, i.e. SkM*, 
Skll, and NLl, are probably not reliable in this region. 
Figure ^ shows that the power of a force for predicting 
total binding energies is fairly independent of its predic- 
tive power for shell effects. Forces with a similar (good) 
description of the smooth trends of binding energies can 
yield rather different magic numbers; compare, e.g., SLy4 
and NL3. 



VII. CONCLUSIONS 

The recent experimental progress in the search for new 
superheavy elements opens a new window for system- 
atic explorations of the limit of nuclear mass and charge. 
Theoretically, predictions in the region of SHE are bound 
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FIG. 8. Macroscopic energy -Emacro.LDM (f^ for the A''=184 isotones as a function of Z calculated for several Skyrme forces 
(lines) and the standard liquid drop model (dots). To illustrate the Z dependence, all energies were normalized to zero at Z=100. 
The bulk parameters of Skyrme forces are given in Table |. The Coulomb-energy constant was assumed to be acoui=0.717 
MeV H in all cases. 



to be extrapolations from the lighter systems. An inter- 
esting and novel feature of SHEs is that the Coulomb 
interaction can no longer be treated as a small pertur- 
bation atop the nuclear mean field; its feedback on the 
nuclear potential is significant. 

The main objective of this study was to perform a de- 
tailed analysis of shell effects in SHE. Since many nu- 
clei from this region are close to the proton drip line, a 
new method of calculating shell corrections, based on the 
Green's function approach, had to be developed. This 
technique was applied to a family of Skyrme interac- 
tions and to several RMF parametrizations. This tool 
turned out to be extremely useful for analyzing the spec- 
tral properties of self-consistent mean fields. 

It has been concluded that both the SHE and RME 
calculations are internally consistent. That is, all 
the Skyrme models employed in this work predict the 
strongest spherical shell effect at iV=184 and Z=124,126. 
On the other hand, all the RME parametrizations yield 
the strongest shell effect at A^=172 and Z=120. It is very 
likely that the main factor contributing to this difference 
is the spin-orbit interaction, or rather its isospin depen- 
dence [BO 59 62-p4] . The role of the spin-orbit potential 
in determining the stability of SHE was posed already 
in the seventies p7| , |88| . The experimental determination 
of the centre of shell stability in the region of SHE will, 
therefore, be of extreme importance for pinning down the 
question of the spin-orbit force. 

Another interesting conclusion of our work is that the 
pseudo-spin symmetry seems to be strongly violated in 
the RME calculations for SHE. As a matter of fact, the 
iV=172 and Z~120 magic gaps predicted in the relativis- 
tic model appear as a direct consequence of pseudo-spin 



breaking. This is quite surprising in light of several recent 
works on the pseudo-spin conservation in RME |^,|9| . 

Einally, from calculated masses we extracted self- 
consistent macroscopic energies. They show a significant 
spread when extrapolating to unknown SHE. This is ex- 
pected to give rise to systematic (smooth) deviations be- 
tween masses and mass differences obtained in various 
self-consistent models. 
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